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Abstract 

We present an algebraic isospin approach towards a more straightforward 
and model- independent determination of the weak phase a (or (j)2) and QCD 
penguin pohution m. B ^ vrvr decays. The world averages of current ex- 
perimental data allow us to impose some useful constraints on the isospin 
parameters of i? — > vrvr transitions. We find that the magnitude of a (or (j)2) 
extracted from the indirect CP violation in vr'^vr" mode is in agreement with 
the standard-model expectation from other indirect measurements, but its 
four-fold discrete ambiguity has to be resolved in the near future. 
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[~T] The major goal of KEK and SLAC 5-meson factories is to test the Kobayashi- 
Maskawa mechanism of CP violation [1] within the standard model and to detect possible 
new sources of CP violation beyond the standard model. An elegant description of CP 
violation in B physics is the unitarity triangle defined by the following orthogonality relation 
of six quark mixing matrix elements in the complex plane [2] : 

VudV:, + KdK; + VtaVl = . (1) 

Three inner angles of this triangle are denoted as a, 13 and 7 by the BaBar Collaboration, 
or equivalently 0i, 02 and 03 by the Belle Collaboration: 



a = 02 = arg 
/3 = 01 = arg 



So far f] has been rather precisely determined from the measurement of CP violation in 
B^j vs B^ J/ipKs transitions [3], and its value (3 ^ 23° is compatible very well with 
the standard-model expectation. The next experimental step is to measure a and 7 to a 
good degree of accuracy at S-meson factories, such that one may cross-check the Kobayashi- 
Maskawa picture of CP violation and probe possible new physics in the S-meson system. 

It is well known that B^ tt+tt", B^ n^n^ and B~ —> tt^tt" decays can be used to 
extract the weak angle a in a model-independent way, because the isospin relation of their 
transition amplitudes allows us to remove the QCD penguin pollution [4]. Note, however, 
that the experimentally-reported branching fractions of S — > tttt decays are all charge- 
averaged: 



+ 2 



B{B', ^ TT+TT-) + B{B', TT+TT-) 

Boo ^ \ [B{B', - + B(^, ^ 



2 

Bo± ^ I [B{B^ + B{B- . (3) 

The world averages of current BaBar [5], Belle [6] and CLEO [7] data on Boo and 

Bo± are listed in Table 1 [8]. In addition, the direct CP- violating asymmetries between 
— > 7r+7r~, — > tt^tt'^, B~ — > 7r°7r~ and their CP-conjugate decays can be defined as 



B{B', TT+TT-) - B{BI TT+TT-) 



B{B^,^n+Tx-)+B{Bl^n+n-) ' 
_ B{B- n^n-) - B{B+ ^ 7r°7r+) 

~ B{B- ^ 77%-) + B{B+ ^ 77%+) ■ ^ ' 
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The world averages of current BaBar [5] and Belle [6] data on C+_, Coo and Ao± are also 
shown in Tabic 1 [8] . These two collaborations have actually measured the time-dependent 



rates of vs 5^ — >• 7r"'"7r decays on the T(45') resonance: 

e-\^t\/ro 

r[5^(At) ^ 7r+7r~] = [1 - S+- sin (Am^^At) + C+_ cos (Am^^At)] , 

rp^(AO ^ = — [1 + Sin (A™,At) - C,. cos (A^.A*)] . (5) 

where Tq is the lifetime of neutral B mesons, and signifies the indirect CP violation 
arising from the interplay between decay and B'^-'b'^^ mixing [9]. A similar time-dependent 
measurement can be done for B'^ vs B^ n^n^ decays, whose rates consist of Coo and iSoo 
corresponding to C+_ and >S+_ in Eq. (5). Only »S+_ has been determined from the KEK 
and SLAC experiments, and its world average [8] is given in Table 1. 

Although the relevant experimental data on direct CP violation remain quite preliminary, 
they can be used to do a quantitative analysis of B nn decays. The main purpose of 
this paper is to recommend an algebraic isospin approach, which allows us to figure out 
the ranges of S — > tttt isospin parameters and to determine the weak phase a (or ^2) from 
»S_|__ and (or) iSoo in a more straightforward and model-independent way. We find that the 
allowed region of a is in agreement with the standard-model expectation from other indirect 
measurements, but its four-fold discrete ambiguity has to be resolved in the near future. 



Under isospin symmetry and in the neglect of electroweak penguin contributions [10], 
the amplitudes of B^ — tt+tt", B^ 7r°7r° and B^ Tr^TT"*" decays (or their CP-conjugate 
processes) form a triangle in the complex plane: 

A(fi° ^ TT+TT") + V2A{B'^ ^ 7r°7r°) = V2A{B+ tt\+) , 

A(S2 ^ TT+TT") + V2A{B'^ ^ TT°TT°) = ^2^(5" ^ TT^TT") . (6) 

The magnitudes of A{B~ — * tt^tt") and ^4(5+ tt°tt+) are identical to each other in 
this safe approximation [4,11], hence the CP- violating asymmetry Ao± vanishes. Table 1 
indicates that the experimental data are in good agreement with the expectation of Ao± ~ 0. 
In terms of the charge-averaged branching fractions in Eq. (3) and the direct CP-violating 
asymmetries in Eq. (4), let us follow Ref. [12] to exphcitly express the parameters 

I I ie _ ^0 
r = |r|e = — , 

A2 

r=\f\e^'^^, (7) 

which stand for the ratios of / = and 1 — 2 isospin amplitudes in B^ — > tt+tt" (or tt'^tt'^) 
and — > TT+TT~ (or tt°tt°) decays. The results are 





(l-FC+_) + 3i3oo 


;i + Coo) -2«Bo± 






(l-C+_)+3Soo 


[l-Coo)-2fi;So± 
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and 



6* = ± arccos 



9 — ± arccos 



6B00 (1 + Coo) - 3B+_ (1 + C+-) - 2kBo± 




(9) 



where k = Tq/t± = 0.921 ± 0.017 [2] denotes the hfetime ratio of neutral and charged B 



\r\ and 9 



would hold, if 



mesons. Eqs. (8) and (9) clearly show that r = r (i.e., |r| 
C+_ = Coo = held. Hence the deviation of r from r is a measure of direct CP violation in 
vs B^^ — > 7r+7r~ and 7r°7r° decays ^. 

The CP-violating parameters S+- and <Soo in Eq. (5) are related to the weak angle a as 
follows [12]: 



'oo 



:i + C+_)Im 



(l + Coo)Im 



p A{B^, 

q A{^,- 



TT+TT") 



= (l + C+_)|i?|sin[2(a; + e)] , 
(l + Coo)|:R|sin[2(o; + e)] , 



(10) 



where q/p^ (^tdytb) I (^tdVtb) denotes the weak phase of B^-B^ mixing [9], and 



2ie 



R= \R\e 



1 -f 

1 — r ' 

2 + r 

2 + r ■ 



(11) 



To be specific, we have 



\R\ 



\R\ 



1 - 2 f cos 61 + 



^1 — 2|r| cos^ + |rp 



4 + 4|r| cos 6* 



4 + 4|r| cos6' 



(12) 



^It is worth mentioning that the difference between ^ = arg[A(B^ — 7r+7r~)/A(B° — 7r°7r°)] and 
99 = aTg[A{B^ TT^TT~) /A{B^ 7r°7r°)] measures the existence of direct CP violation too [12]. 
With the help of Eqs. (3), (4) and (6), one may obtain 

2kSo± - e+- (1 - - 2^00 (1 - Coo) 



COS<^ : 

COS if ■ 



2v/2fi+-^?oo (1 - C+_) (1 - Coo) 
2Kfio± - (1 + C+_) - 2^00 (1 + Coo) 



2v/2e+_^oo(l + C+_)(l + Coo) 

Obviously, cos ip = cos (j) would hold if both and Coo were vanishing. Note that the notations 

of direct CP-violating asymmetries in Ref. [12] are A-\ = — C+_ and .4oo = — Coo- 
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and 



1 

G = + - arctan 

© = — arctan 
2 



r sin Q — 


r sin 6' — r 


r 


sin(^ — 






1 — |r| cos^ 


— |f| cos^ + 


r| 


r| cos(^ 






2 r| sin^ - 


- 2|f| sin^ + 


|r| 


|f| sin(0 


-^) 1 



4 + 2|r| cos6' + 2|f| 008 6* + |r||r| cos(6' - Q) 



(13) 



If there were no direct CP violation in 5 — > tttt transitions (i.e., C+_ = Cqo = or f = r), we 
would arrive R — R — \ from Eqs. (11)-(13). In this case, Eq. (10) would be simphfied 
to »S+_ = »Soo = sin2Q;. It is therefore necessary to pin down C+_ and Cqo to a reasonable 
degree of accuracy, in order to extract the weak angle a from »S_|__ and (or) »Soo. 



Now we carry out a numerical analysis of S — > tttt decays and CP violation by using 
the isospin formulas obtained above and current experimental data listed in Table 1. The 
Icr, 2(7 and Scr confidence regions of (|r|, |f|), (cos^,cos^), (|i?|, |i?|), (9,9) and (q;,iS+_) 
are shown in Figs. 1-3, while the central (best-fit) values of these parameters are given in 
Table 2. Some discussions are in order. 

(1) The moduli |r|, |f|, \R\ and \R\ can be determined without any discrete ambiguity 
(see Fig. 1 for illustration). The difference between |r| and |f| is quite obvious, although it 
remains possible for |r| = |r| to hold at the 3(T level. In contrast, it is likely to have \R\ = \R\ 
at the la level, but their best-fit values arc different from each other. Direct CP violation 
is therefore expected to manifest itself in B nn transitions. 

(2) Although cos 9 and cos 9 can be uniquely determined, 9 or 9 involves two-fold ambi- 
guity. The difference between 9 and 9 is apparent, but cos^ = cos^ remains possible at the 
1(7 level. More precise experimental data will allow us to fix r and f, both their moduli and 
their phases, to a better degree of accuracy. 

(3) The two- fold ambiguity of 6' or 6' leads to the four-fold ambiguity of 9 or 9, as 
illustrated in Fig. 2 and Table 2. An interesting feature of our results is that the signs of 
9 and 9 are essentially opposite. Although it is possible to have 9 = 9 at the 2a level 
for case (A) or (D), we find that 9 7^ 9 holds even at the 3(7 level for case (B) or (C). 
In particular, the possibility of 9 = and (or) 9 = is strongly disfavored, implying the 
presence of QCD penguin pollution or direct CP violation in i? — > tttt decay modes. The 
typical values of 9 and 9 are 9 ~ ±19° and 9 ~ =f26°, as shown in Table 2. 

(4) Because of the four- fold ambiguity associated with 9 or 9, the result of a determined 
from iS+„ involves the four-fold discrete ambiguity too, as illustrated in Fig. 3. Table 2 tells 
us that the central values of a can be a ~ 122° (A), 135° (B), 86° (C) or 95° (D). This result 
is certainly in agreement with the standard-model expectation of a (i.e., a ~ 90°) from other 
indirect measurements [2]. Once the CP- violating asymmetry Sqq is also measured, it will 
be possible to completely or partly remove the discrete ambiguity of a [4]. Then one may 
constrain the weak phase a and QCD penguin pollution in S — > tttt decays at a much better 
confidence level. 

It is worth remarking that the validity of our isospin analysis relies on the assumption 
of negligible electroweak penguin effects. The electroweak penguin contribution to S — > tttt 
decays is in general expected to be insignificant [10]. This expectation would be problematic 
or incorrect, if Aq± ^ were experimentally established [12]. Note also that final-state 



5 



interactions in S — > tttt transitions consist of both elastic tttt tttt rescattering and some 
possible inelastic rescattering effects ^. Whether the latter is negligibly small or not remains 
an open question [15]. To answer this question requires more precise measurements of both 
branching fractions and CP-violating asymmetries of 5 —> tttt decays. 



We have presented an algebraic isospin analysis of rare B tttt decays by taking 
account of the fact that the experimentally-reported branching fractions are charge-averaged 
and large direct CP violation may exist in them. This approach is more straightforward 
than the originally-proposed geometric approach, from which the weak phase a (or ^2) and 
QCD penguin pollution are determined through the reconstruction of two isospin triangles. 
Therefore, our method is expected to be very useful to analyze the future experimental data 
on B ^ TTTT transitions and CP violation in a model-independent way. 

Although the present experimental data (in particular, those on direct and indirect CP 
violation in S — > tttt decays) are not sufficiently precise, they can impose some instructive 
constraints on the parameter space of QCD penguin effects. Furthermore, we find that the 
allowed region of a (or is actually in agreement with the standard-model expectation 
from other indirect measurements. To resolve the four-fold discrete ambiguity associated 
with the magnitude of a (or ^2) determined from the indirect CP-violating asymmetry 
in TT'^TT^' mode, a measurement of the similar CP-violating asymmetry in 7r°7r° mode is 
necessary. We expect that more accurate measurements of such charmless B decays will 
help us to test the consistency of the Kobayashi-Maskawa mechanism of CP violation ^ and 
to probe possible new physics beyond the standard model. 

This work was supported in part by the National Nature Science Foundation of China. 



^The transition amplitudes of B'^ vr'^vr", B'^ tt^tt^ and B^ tt^tt^ decay modes (or their 
CP-conjugate processes) may still form an isospin triangle in the complex plane, even if the inelastic 
TTTT # DD rescattering effects are taken into account [13]. In this complicated case, however, a 
model-independent determination of a from <S_| and <Soo would be rather difficult [14] . 

^For example, the relationship sin a/ sin /3 = \Vcd/Vud\/\Vub/Vcb\ [16], which holds as a straight- 
forward result of the unitarity-triangle defined in Eq. (1), can be numerically tested with more 
accurate data of a and |V^ub/Kb|- 
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TABLES 



TABLE I. The world averages of current experimental data on the charge-averaged branch- 
ing fractions (H^ , ^Boo^ Bo±), direct CP-violating asymmetries (C^ , Cqq, Ao±) and indirect 

CP-violating asymmetries 5oo) of B —>■ tttt decays [8]. 





World average 


B+- 


(4.6 ± 0.4) X 10-6 


Boo 


(1.51 ±0.28) X 10-6 


l3o± 


(5.5 ±0.6) X 10-6 


C+- 


-0.37 ±0.11 


Coo 


-0.28 ± 0.39 


^o± 


-0.02 ± 0.07 


S+- 


-0.61 ±0.14 


«5oo 





TABLE IL The central values of eight isospin parameters (|r|, |f |; 9, 6; \R\, \R\; 6, 6) and the 
weak phase a constrained from the world averages of current BaBar and Belle data on B ^ tttt 
decays [8], where we have taken into account the two- fold ambiguity associated with and 6 as 
well as the four-fold ambiguity associated with 0, 6 and a. 





Case (A) 


Case (B) 


Case (C) 


Case (D) 


\r\ 


0.6 


0.6 


0.6 


0.6 


\r\ 


1.7 


1.7 


1.7 


1.7 


9 


±180° 


-180° 


±180° 


-180° 


9 


+120° 


±120° 


-120° 


-120° 


\R\ 


1.5 


1.5 


1.5 


1.5 


\R\ 


1.3 


1.3 


1.3 


1.3 


Q 


-18° 


-20° 


±20° 


±18° 


Q 


±25° 


±27° 


-27° 


-25° 


a 


122° 


135° 


86° 


95° 
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FIGURES 
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\R\ 

FIG. 1. The la, 2a and 3a confidence regions of (|r|, |r|), (cos 6, cos 9) and (|i?|, |i?|) parameters, 
constrained by the isospin relations and current experimental data. 
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FIG. 3. The la, 2a and 3a confidence regions of a with four-fold discrete ambiguity, extracted 
from current experimental data on S-\ 
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